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Abstract. A search for stellar clusters has been carried out in 18 nearby spiral galaxies, using archive images 
from the Wide Field Planetary Camera 2 on board the Hubble Space Telescope. All of the galaxies have previously 
been imaged from the ground in UBVI. A catalogue of structural parameters, photometry and comments based 
1 on visual inspection of the clusters is compiled and used to investigate correlations between cluster structure, 

environment, age and mass. Least-squares fits to the data suggest correlations between both the full-width at 
half- maximum (FWHM) and half-light radius (R B a) of the clusters and their masses (Af) at about the 3<r level. 
Although both relations show a large scatter, the fits have substantially shallower slopes than for a constant- 
density relation (size oc M 1//3 ). However, many of the youngest clusters have extended halos which make the 
-Reff determinations uncertain. There is no evidence for galaxy-to-galaxy variations in the mean cluster sizes. In 
particular, the mean sizes do not appear to depend on the host galaxy star formation rate surface density. Many 
of the youngest objects (age < 10 7 years) are located in strongly crowded regions, and about 1/3-1/2 of them are 
double or multiple sources. The HST images are also used to check the nature of cluster candidates identified in 
' a previous ground-based survey. The contamination rate in the ground-based sample is generally less than about 

20%, but some cluster identifications remain ambiguous because of crowding even with HST imaging, especially 
for the youngest objects. 
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^ ' 1. Introduction ber of unusually (by Milky Way standards, at least) bright 

and massive young clusters (Shapley & Nail 1951. Hodge 
In previous papers (Larsen & Richtler EH3 Larsen [M3 1961 Richtler gggg christian & Schommer EH1 [1988; 
£ ■ Larsen & Richtler 2000. hereafter Paper I, I and III), we chandar ^ &1 ^ ^ paper ffl wg conduded that the 

rTN , have studied populations of young stellar clusters m the , . . r , , . , in 1>JT , 

5— I . . , . . . , . . mam driving factor behind these differences seems to be 

rrt disks of nearby spiral galaxies using ground- based imag- ,, , , ,. , /0 t^ti\ r , L i ± i n i 

. . . . i • iii- i the star formation rate (Si R) of the host galaxy. Galaxies 

me. I hat work was motivated by a desire to understand .,, , . , OT ™ / . , n ,, f 

° ,ii with high Si Rs (per unit disk area) apparently lorm a 

why some galaxies host young stellar clusters which are , - , . ,,, . . . . , , , , 

J ° . larger fraction of their stars m massive, bound clusters, 

significantly more luminous (and, presumably, more mas- mi. ru- ii i i j. i 

° . J K ' f J \ I he presence of highly luminous clusters m galaxies with 

sive) than open clusters m the Milky Way. Well-known , . , , , , , , . c , 

, r , , . high Si Rs may be — at least partially — a size-of-samplc 

examples of galaxies with rich populations of luminous *r j. j j. it. • i. i i w • u i 

r ° . , r ,. , effect, due to the rich cluster populations m such galaxies 

young clusters include a number of merger galaxies and ^ &1 ^ Larsenr2lm 

starbursts (see e.g. compilation m Whitmore 2003), but 

there are also some relatively "normal" galaxies such as The study of stellar clusters is intimately linked to that 

the Large Magellanic Cloud and M33 which host a num- of star formation in general. Observations show that a 

large fraction, if not the vast majority, of all stars are born 

* Tables El H |U and □ are only available in clusters (e.g. Carpenter I21TOI Lada fc Ladal5003). This 

in electronic form at the CDS via anonymous , -lii^n i_ jj j i 

, „ ,„„ . does not, however, imply that all embedded clusters are 

ftp to cdsarc.u-strasbg.fr (130.79.125.5) or via . . ...... 

rrrr — m r — — r — j. , . , ^ ; . , — n dynamically bound entities which survive emergence from 

http://cdsweb.u-strasbg.fr/Abstract.ntml J J ° 

** Based on observations obtained with the NASA/ESA their native molecular cloud cores and become observable 

Hubble Space Telescope, obtained at the Space Telescope at optical wavelengths. Lada & Lada H2003J estimated 

Science Institute, which is operated by the Association of that less than 4-7% of local embedded clusters survive to 

Universities for Research in Astronomy, Inc., under NASA con- become bound clusters of Pleiades age (~ 10 s years), but 

tract NAS 5-26555 this number may depend on environment. In Paper III we 



2 



S. S. Larsen: Young Stellar Clusters in Spirals 



noted a steady increase with host galaxy area-normalised 
star formation rate (Esfr) in the fraction of [/-band light 
originating from clusters, ranging from well below 1% in 
galaxies with very low SFRs (like IC 1613) to several per- 
cent in starbursts. Meurer et al. 1^)95) found that on av- 
erage about 20% of the UV light in a sample of starburst 
galaxies comes from young clusters. De Grijs et al. (2003) 
estimated that as much as ~ 70% of the B-band light 
in the tidal tails of the "Tadpole" and "Mice" galaxies 
originates from young clusters or compact star-forming re- 
gions. In most "normal" star-forming galaxies, the young 
clusters contribute about 1% of the £/-band light. Whether 
or not an embedded cluster remains gravitationally bound 
depends on the star formation efficiency within the proto- 
cluster cloud, as well as the timescale on which the gas 
is expelled (Elmegreen 1983 Kroupa 2001). It may there- 
fore be more appropriate to view the fraction of optically 
visible young stars associated with clusters as a survival 
frequency, associated with the star formation efficiency, 
than a cluster formation efficiency per se (which is prob- 
ably always close to 100%). 

While care must be taken when interpreting the above 
results, due to possible differences in the age distributions 
of the clusters and/or field stars, completeness limits, etc., 
it seems clear that the Solar neighbourhood samples only 
a small part of the conditions under which star forma- 
tion takes place in the Universe. Fortunately, there are 
several star- forming galaxies available within a few Mpc, 
spanning a range in SFRs, morphological type etc., which 
can be studied in considerable detail with a combination 
of ground-based and space-based techniques. The origi- 
nal sample of 21 nearby spirals analysed in Papers I— III 
has since been augmented by an additional handful of 
galaxies observed with the 3-m Shane telescope at Lick 
Observatory (see Larsen 2002). However, on the ground- 
based images, clusters were only marginally resolved, and 
although significant efforts were made to weed out the 
most obvious contaminants, the cluster lists in Paper II 
should only be taken as provisionary. Many of the galax- 
ies have now been imaged with the Wide Field Planetary 
Camera 2 (WFPC2) on board the Hubble Space Telescope 
(HST) for a variety of reasons, often with multiple point- 
ings. The WFPC2 imaging not only provides a welcome 
check of the true nature of the sources identified as cluster 
candidates from the ground, but also allows the structure 
and immediate environment of individual clusters to be 
examined in much greater detail. 

Several studies have indicated a remarkable unifor- 
mity in the sizes of stellar clusters over a wide range of 
masses, environments, ages and metallicities. The most 
robust measure of cluster size is the half-light or "effec- 
tive" radius (R c s) which is expected to remain relatively 
stable over the lifetime of a cluster (Spitzer 1987). For 
Galactic globular clusters, i? e ff and luminosity are uncor- 
rected, although there is a trend of increasing cluster size 
with galactocentric distance (Rqc', van den Bergh et al. 
I1991J) . Using data from the compilation by Harris (1996), 
the median i? c ff is 3.0 pc. Similarly, the diameters and 



masses of Galactic open clusters show no strong correla- 
tion, with typical sizes only slightly smaller than those of 
globular clusters (Janes et al. 1988). For young clusters 
in the "Antennae" merger, Whitmore et al. (1999) found 
mean effective radii of 4±1 pc. Zepf et al. ill 999) estimated 
half-light radii of 5-10 pc for clusters in NGC 3256, per- 
haps slightly larger than for the Antennae, but again with- 
out any strong size-luminosity correlation. For globular 
clusters around early-type galaxies, typical effective radii 
are again 3-4 pc with no clear size-mass correlation (e.g. 
Kundu & Whitmore 2001). The lack of a significant mass- 
size relation is puzzling, since one might a priori expect 
a cluster to form once the parent gas cloud reaches a cer- 
tain density, independent of the total mass. If this initial 
gas density is reflected in the stellar density of the result- 
ing cluster, one might naively expect the radius to scale 
with mass (M) roughly as i? e ff oc M 1 / 3 . However, this 
is not what has generally been reported. From the above 
examples, it appears that star clusters typically have ef- 
fective radii of a few (~ 3) pc, with a scatter of perhaps 
1-2 pc. Exceptions are found, however, including the faint 
"Palomar" -type globular clusters in the outer part of the 
Galactic halo, and the "faint fuzzy" clusters recently dis- 
covered in a couple of nearby SO-type galaxies (Larsen & 
BrodieEOOni Brodie & Larsen l2T)Tl2"|l . which have larger 
effective radii (^ 10 pc). 

In this paper, the cluster candidates identified from 
the ground are first re-examined on archive WFPC2 im- 
ages. Additional cluster candidates are then identified on 
the WFPC2 images and combined with ground-based pho- 
tometry to produce a catalog of structural parameters and 
photometry for a sample of clusters. Relying on ground- 
based photometry limits the sample to relatively bright 
objects, but has the advantage of providing uniform pho- 
tometric coverage of all clusters (even if crowding effects 
are more severe than in the HST data). In particular, most 
of the HST datasets do not include imaging in a U -band 
equivalent filter, which is essential for age-dating the clus- 
ters. However, the HST photometry may still be useful 
for some purposes and aperture photometry in an r = 0'.'5 
aperture is presented for the available bandpasses in a 
separate table for each cluster candidate. Each entry in 
the catalogue also contains various notes on the degree of 
crowding, close neighbours etc. based on a visual inspec- 
tion of the candidates. 

Because the HST images cover limited sections of the 
galaxies and span a huge range in exposure times and 
filters, the cluster sample presented here still cannot be 
considered complete in any sense. The completeness is a 
complicated function of crowding, cluster size, underly- 
ing surface brightness, exposure time in the HST images, 
bandpass, seeing in the ground-based data, galaxy dis- 
tance, and probably many other factors which would be 
next to impossible to model in a satisfactory way. It should 
also be emphasized that what is presented is still a list of 
cluster candidates, which might contain contaminants (e.g. 
background galaxies). A definitive list of bona-fide clusters 
would require spectroscopic follow-up, but such an effort is 
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Table 1. Exposures. A(a,6) indicates the offsets in right 
ascension and declination between the coordinate systems 
of the WFPC2 frames and the ground-based data (tied to 
the USNO catalogue). The offsets are given in sec/15 for 
right ascension and in arcseconds for declination (i.e. the 
cos 5 factor has not been applied to convert the offsets in 
a to true arcsecs). 



Table 1. Exposures - continued 
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2. Data and Reduction 

The search for HST archive data was concluded in October 
2002 and only includes WFPC2 datasets which had been 
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beyond the scope of this paper and left for future studies. 
With this in mind, it is hoped that the catalogue may still 
provide a useful basis for further studies. As an example, 
it is used in Section 0] to investigate trends with age and 
mass in the cluster sizes, degree of crowding, and shape 
parameters. 



publicly released up until that time (Table HJ. ACS data 
were available for a few galaxies, but have been excluded 
in order to allow a relatively simple and homogeneous re- 
duction procedure. 

Given the large volume of data, a fairly high degree 
of automatization had to be incorporated in the reduc- 
tion procedures. When several exposures were available 
for a given field and filter, these were combined with the 
CRREJ task in the STSDAS package in IRAF 1 . In most 
cases, no shifts in the image coordinate systems were re- 
quired before combination but when necessary, such shifts 
were applied using the IMSHIFT command in IRAF. For 
each WFPC2 pointing and each filter, objects were then 
detected with the DAOFIND task in DAOPHOT (Stetson 
1987) running within IRAF. Objects with saturated pix- 



1 IRAF is distributed by the National Optical Astronomical 
Observatories, which are operated by the Association of 



Universities for Research in Astronomy, Inc. 
with the National Science Foundation 



under contract 
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els within a radius of 5 pixels were rejected. The detected 
sources were fitted with the ISHAPE profile-fitting algo- 
rithm (Paper II). ISHAPE models each source assuming 
an analytic model for the intrinsic profile of the source, 
convolved with the HST/WFPC2 point spread function 
(PSF). The FWHM of the analytic model is iteratively 
adjusted until the best fit to the data is obtained. The 
initial round of profile fitting was carried out with a fit- 
ting radius of r = 5 pixels and assuming a model of the 
form 

P(r) oc (l + (r/r c ) 2 p (1) 

with a = 1.5. Profiles of this type were shown by Elson 
et al. I|1987[l to provide good fits to young LMC clusters, 
and are hereafter referred to as "EFF" models (but note 
that the exponent a in Eq. Q corresponds to 7/2 in the 
definition by Elson et al. 1987). The WFPC2 PSF was 
modelled using version 5 of the TINYTIM package (Krist 
& Hook 1997), including a convolution with the "diffusion 
kernel" to account for pixel-to-pixel charge diffusion. PSFs 
were generated for each filter, and automatically selected 
by the reduction scripts to match the bandpasses used for 
the observations. 

Because the main aim was to study the structure of 
spatially resolved objects, only objects which were de- 
tected at > lOcr above the background noise were included 
in the analysis. Fainter objects would generally have too 
low S/N for reliable size measurements, and would have 
increased the already substantial computing time required 
to fit the spatial profiles. Because of the vastly different ex- 
posure times, different bandpasses, background levels etc., 
the 10-er detection threshold does not translate into a well- 
defined completeness limit in terms of magnitude. A total 
of 82000 sources were detected and fitted, requiring a few 
days of CPU time on a 1.5 GHz Pentium PC. At this stage, 
many objects appeared several times in the source list, be- 
ing included in multiple HST pointings and/or bands. 

After the initial round of profile fitting, the HST object 
lists were matched with the photometry data files from 
the ground-based surveys (all details concerning the re- 
duction of the ground-based data are given in Paper II). 
Not only cluster candidates previously identified as such 
in the ground-based surveys, but all point-like sources in 
each ground-based CCD frame for which photometry was 
available, were matched. The matching was done by con- 
verting pixel coordinates measured on the WFPC2 im- 
ages to celestial coordinates, using the METRIC task in 
the STSDAS package in IRAF. These coordinates were 
then matched with the coordinates of objects measured 
on the ground-based CCD images, tied to the US Naval 
Observatory meridian circle catalogue as provided by the 
ESO SkyCat Tool (Monet et al. JHffito . In many cases, 
there were clear systematic offsets between the WFPC2 
and USNO coordinate systems. These offsets were deter- 
mined by displaying each individual WFPC2 frame and 
then marking the ground-selected cluster candidates con- 
tained within that frame. Because these would typically 
be among the brightest objects in the WFPC2 frames, 



Table 2. Explanation of comment codes 

a double/multiple object of comparable brightness 

b fainter neighbours within tf.'5 

c companions of similar brightness between 0'.'5 < r < l'.'5 

d elongated 

e chain of compact sources 

f irregular 

g crowded 

h in group/association of more than 5 bright stars/clusters 

i irregular background 

j saturated star 

k near edge of CCD 

1 star? 

m galaxy? 
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12 3 4 5 

FWHM (pixels) 

Fig. 1. Histogram of FWHM values for objects fitted dur- 
ing the first round of profile fitting. The narrow sequence 
of unresolved objects with FWHMwO is clearly seen. The 
vertical dashed line indicates the criteria used to select 
cluster candidates. 

and the offsets in general relatively minor (<J 1"), iden- 
tifying the cluster candidates in the WFPC2 frames was 
usually unproblematic. The offsets between the WFPC2 
and USNO coordinate systems are listed in Table ^ f° r 
each dataset. Note that large offsets were found for the 
exposures belonging to snapshot programme 5446. This 
is probably because these datasets were guided using the 
"Gyro Hold" mode, which provides less accurate pointing 
and tracking than the Fine Guidance Sensors on HST. A 
few WFPC2 pointings had no ground-based cluster candi- 
dates and were excluded from further analysis. For those 
pointings, no offsets are given in Tabled 

The initial round of profile fitting served to separate 
point-like sources (likely stars) from extended ones (clus- 
ter candidates). Fig. ^ shows the distribution of intrinsic 
FWHM values (in pixels). There is a very narrow peak of 
objects with FWHM values close to 0, although there is 
no sharp division between resolved and unresolved objects. 
For further analysis, a total of 3100 sources with S/N>50 
(measured within the 5 pixels aperture on the WFPC2 
images), FWHM>0.2 pixels (dashed line in Fig. and 
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a (double/multiple): 

nl313-540, n4258-1132, n5194-477 



b (fainter neighbours): 
nl313-304, n4258-1051, n5194-239 



-V 



c (companions within l'.'5): 
nl313-471, n4258-454, n5194-121 



d (elongated): 
n2835-625, n3184-277, n5194-351 




■ ' ■ 



e (chain): 

n5194-424, n5194-906, n628-1895 



f (irregular): 

nl313-279, n2835-890, n7424-142 



g (crowded): 

n5194-897, nl313-508, n2997-558 

■ , - -* ■ 

h (in group of more than 5 stars/clusters): 
n4258-776, n5194-433, n6946-1489 



. *- 



i (irregular background): 
n4258-1060, n4258-794, n5194-394 







Jt • « 








■ • ■ . 




Fig. 2. Illustrations of comment codes a-i. Please see 
Table |2 for the full comments associated with each code. 
Each panel shows a 4" x 4" section around the cluster 
candidate. 



Fig. 3. Examples of objects classified as Type 2 (uncer- 
tain) 



N3621-513: 

F336W F439W F555W F814W 



N3621-620: 

F336W F439W F555W F814W 




Fig. 4. Two examples of how the appearance of cluster 
candidates can change dramatically depending on band- 
pass. 



ground-based photometry were selected. For comparison, 
the (undersampled) PSF of WFPC2 itself has a FWHM of 
about 1.5 pixels. At a typical distance of 5 Mpc, the size 
cut corresponds to a FWHM of 0.5 pc or a core radius 
of about 0.25 pc. Many of these sources still represented 
multiple observations of the same object. Objects satisfy- 
ing these criteria were again fitted, but this time allowing 
the power-law index (a) of the EFF profiles to vary as a 
free parameter. In order to better constrain a, a fitting 
radius of 10 instead of 5 pixels was chosen for this second 
round. The choice of a 10 pixels fitting radius represents 
a compromise between reasonably accurate constraints on 
the fit parameters, reducing the effects of a non-uniform 
background, and computing time. Each source was fitted 
three times, using different initial guesses for a (olq — 1.0, 
1.5 and 2.0). The values 1.0 and 1.5 roughly bracket those 
typically found for real star clusters (e.g. Elson et al.[l987 
and Section l4.3|) . while an additional more extreme value 
of ao = 2 was included in order to reduce any a-priori bias 
in the measurements, in case some clusters have steeper 
slopes (although tests of the ISHAPE algorithm suggest 
that any dependence of the fitted parameters on the input 
guesses is relatively minor; see Sec 12.21 below'). 

Finally, all cluster candidates fitted in the second 
round were visually inspected. In addition, all cluster can- 
didates identified in the ground-based surveys and with 
HST imaging were inspected regardless of the S/N in 
the HST images. The inspection was done by display- 
ing a 4" x 4" section of all images of each cluster can- 
didate in an IDL widget, which also contained a num- 
ber of check-boxes corresponding to the comment codes 
given in Table Each object was also assigned to one 
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Table 4. Additional comments (only 5 sample entries 
given) 

ID Comments 

nl313-356 identification uncertain 

nl313-409 tracking error 

nl313-411 cluster members resolved? 

n2835-677 very elongated halo 

n2997-616 bad pixels 



of three types: 1) likely/certain cluster, 2) uncertain clas- 
sification, 3) likely /certain non-cluster. A flag was set if 
the ISHAPE fit was unlikely to have resulted in mean- 
ingful structural parameters, as judged from the degree of 
crowding etc. The inspection process was repeated three 
times, each time displaying the cluster candidates in a dif- 
ferent (random) order. At the end, the final object type 
and ISHAPE "Fit-OK" flag were determined as the most 
pessimistic of the three individual estimates. A comment 
code was included in the final list if selected in at least 2 
out of the 3 inspection rounds. 

Clearly, by its very nature any visual inspection in- 
volves an element of subjectivity. Very often, more than 
one comment applied to a given object. In such cases, all 
applicable comments were selected and included in the 
catalogue. Furthermore, resolution (and thus distance) ef- 
fects may change the visual appearance of certain objects 
and thereby e.g. cause an object classified as "double" in a 
nearby galaxy to be labelled "elongated" in a slightly more 
distant one. In order to illustrate roughly what the indi- 
vidual comment codes in Table represent, Fig. shows 
4" x 4" sections around some clusters that are typical for 
each comment code. Fig.|21shows a few examples of objects 
classified as "uncertain" . Another difficulty is that objects 
can appear quite different in different bandpasses. As an 
example, Fig. 01 shows two objects observed in F336W, 
F439W, F555W and F814W. The F336W and F8I4W im- 
ages of N362 1-620 are hardly recognisable as the same ob- 
ject — the F336W image just shows what might be a loose 
association of stars, while the F814W image shows two 
fairly well-defined, compact sources. Conversely, N3621- 
513 looks compact and symmetric in F336W, while two 
nearby neighbours (perhaps red giants or supergiant stars) 
appear in F814W. If multi-wavelength data were available 
for all objects, one might adopt a consistent strategy for 
dealing with this problem, but in many cases data were 
only available in one band. It should therefore be borne in 
mind that the visual comment codes (and even the profile 
fits) depend on the bandpass in a way that is not easily 
predictable. 

The final list of cluster candidates includes 1358 ob- 
jects (TableOJ). Columns (l)-(8) contain information from 
the ground-based data: Coordinates, UBVI photometry 
and projected galactocentric distance (in arcmin). The 
photometry has been corrected for Galactic foreground 
extinction using the Schlegel et al. (1998 ) values and the 
extinction law in Cardelli et al. 1)198911 . Note that the red- 



dening corrections differ from those in Paper I— III, where 
Burstein & Heiles l|1984fl values were used. In Paper I— 
III we used a relatively large aperture radius of 8 pix- 
els (T/5 and 3'.'2 for the NOT and Danish 1.54 m data, 
respectively) to avoid possible systematic effects in the 
integrated magnitudes due to the extendedness of the ob- 
jects. However, it is now clear that most clusters are com- 
pact enough that this is not a major source of concern in 
ground-based imaging, and in the present paper I therefore 
use a smaller aperture radius of 4 pixels, the same as for 
the colours, for the ground-based magnitudes. Information 
derived from HST images is listed in columns (9)-(16): 
The number of individual detections of each cluster (where 
one "detection" is defined as the presence of the cluster in 
an image taken through a given filter under a given pro- 
gramme), N, is in col. (9), followed by the full- width at 
half maximum (FWHM) of the cluster profile derived from 
the ISHAPE fits, the exponent a, the effective (half-light) 
radius R c s, the minor/major axis ratio (x/y), the object 
type (Col. 14), Fit-OK flag (Col. 15) and comments. 

Unlike the classical King profiles fKing Hl)f)2l IT966). 
the EFF models have no finite radius, and for a < 1 
the volume contained under the profile is infinite. For a 
only slightly larger than unity, the total volume converges 
very slowly, resulting in unrealistically large R c g values. 
To cope with these difficulties, the i? c ff values in Table |21 
are computed for a finite outer radius of 50 pc, beyond 
which the luminosity profiles of young clusters become 
difficult to trace even in nearby galaxies such as the LMC 
(e.g. Elson et al. I1987[) . However, it is important to note 
that the estimates of R e g are generally based on an ex- 
trapolation of the luminosity profiles beyond the fitting 
radius, and carry significant uncertainties especially when 
a 5a 1. Instead of listing the FWHM, I could have given 
the core radius r c , since both are always defined. The rea- 
son for listing FWHM is that there is some ambiguity 
in the definition of the core radius. Some authors define 
it as FWHM/2, but it may also be defined e.g. as the 
scale radius r c in Eq. Jy. In order to avoid confusion, 
I will simply use the term FWHM rather than "core ra- 
dius" throughout the remainder of this paper. When more 
than one exposure was available for a cluster candidate, 
the shape parameters in Table|21were obtained by weight- 
ing the measurements on each exposure by its S/N. As 
discussed above, the morphology of clusters can be quite 
wavelength-dependent, but cases where the determination 
of shape parameters is particularly uncertain can gener- 
ally be recognized by the 'Fit-OK' flag in Column (15) of 
Table El 

Errors were estimated as follows: For each exposure, 
the error on the shape parameters derived from that expo- 
sure were estimated as the standard deviation of the three 
individual fits. If only one exposure was available, this is 
the error listed in Table |21 When several exposures were 
available, the errors in Tableware the estimated standard 
errors on the mean of the weighted average. Some addi- 
tional comments for a few objects, which did not fit into 
the codes in Table are listed in Tabled 
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Table 3. All clusters with S/N>50 on HST images. Photometry is from ground-based data while morphological data 
are from HST/WFPC2 images. The hrst 10 rows of the table are reproduced here; the full table (1358 rows) is only 
available in electronic form. 



ID 


R.A. 


(2000.0) 


Dccl. (2000.0) 


V 


U — B 


B — V 


V — I 


R GC 








JV 


FWHM (pc) 


Oi 


Rcff (pc) 


x/y 


Code 


Fit ok? 


Comments 


(1) 




(2) 


(3) 


(0 


(5) 


(6) 


(7) 


(8) 








0) 


(10) 


(n) 


(12) 


(13) 


(14) 


(15) 


(16) 


nll56-486 


2:59:41.29 


25:14:15.34 


16.97 ± 0.01 


— 0.96 ± 0.01 


0.44 ± 0.01 


— 0.53 ± 0.02 


0.30' 








1 


8.55 ± 0.81 


0.95 ± 0.08 


12.42 it 1.24 


0.63 it 0.04 


1 


N 


a,f 


nll56-583 


2:59:42.85 


25:14:28.68 


17.41 ± 0.01 


—0.29 ± 0.01 


0.29 ± 0.01 


0.65 ± 0.01 


0.23' 








1 


2.03 it 0.55 


0.88 ± 0.05 


10.59 ± 1.11 


0.97 it 0.03 


1 


Y 




nl313-180 


3: 


18:15.12 


— 66:31:55.29 


20.67 ± 0.04 


0.11 ± 0.15 


0.67 ± 0.07 


0.98 ± 0.06 


2.07' 








1 


5.35 ± 0.11 


0.89 ± 0.01 


13.57 ± 0.53 


0.95 it 0.03 


1 


Y 




nl313-199 


3: 


18:17.49 


— 66:31:36.73 


20.31 ± 0.03 


0.03 it 0.08 


0.18 it 0.04 


0.49 ± 0.05 


1.77' 








1 




o no 4- n k& 
*z . uy it u.oo 


o.y/ it u . Dy 


U. (o it U . Uo 








nl313-228 


3: 


18:18.46 


-66:31:22.76 


19.30 ± 0.02 


0.23 ± 0.06 


0.42 ± 0.03 


0.57 ± 0.02 


1.56' 








1 


1.15 ± 0.47 


5.64 ± 5.24 


0.52 ± 0.16 


0.53 it 0.39 


1 


Y 




nl313-234 


3: 


18:07.67 


-66:31:19.06 


17.87 ± 0.01 


-0.82 ± 0.01 


0.34 ± 0.01 


-0.78 it 0.02 


1.66' 








1 


9.24 ± 2.37 


0.87 ± 0.60 


19.43 ± 6.90 


0.92 ± 0.02 


2 


N 




nl313-239 


3: 


18:17.59 


-66:31:12.41 


19.86 ± 0.04 


-0.07 ± 0.08 


0.26 it 0.06 


0.38 it 0.08 


1.37' 








4 


6.34 ± 0.46 


1.29 ± 0.32 


8.90 ± 1.09 


0.90 ± 0.02 


1 


Y 




nl313-248 


3: 


18:16.05 


-66:31:06.91 


18.85 ± 0.03 


1.21 ± 0.24 


1.03 it 0.07 


1.11 ± 0.04 


1.27' 








1 


2.14 ± 0.22 


70.80 ± 90.78 


0.67 ± 0.05 


0.22 ± 0.03 


2 


Y 




nl313-249 


3: 


18:12.45 


-66:31:06.38 


19.95 ± 0.06 


0.10 ± 0.10 


0.20 it 0.07 


0.60 ± 0.08 


1.29' 








3 


5.18 ± 0.26 


1.25 ± 0.08 


7.09 ± 0.70 


0.92 ± 0.03 


1 


Y 




nl313-259 


3: 


18:15.55 


-66:31:02.25 


20.71 ± 0.14 


-0.31 ± 0.19 


0.26 it 0.19 


0.71 ± 0.18 


1.19' 








2 


1.40 ± 0.42 


1.41 ± 0.20 


1.44 ± 0.05 


0.53 it 0.03 


1 


Y 





Table 5. HST photometry for the same clusters listed in Table|3J All magnitudes are in the STMAG system, measured 
in a C/5 aperture and applying an aperture correction of —0.2 mag. No corrections for foreground reddening have been 
applied. Only the hrst 10 rows are reproduced here. 



ID 



F300W 



F336W F439W 



F450W 



F547M F555W 



F606W 



F814W 



nl 156-486 
nl 156-583 
nl313-180 
nl313-199 
nl313-228 
nl313-234 
nl313-239 
nl313-248 
nl313-249 
nl313-259 



18.29 ±0.01 
19.68 ± 0.04 



19.51 ±0.01 

20.38 ± 0.02 
19.32 ± 0.01 
20.28 ± 0.02 



21.11 ±0.01 
21.05 ± 0.01 
19.66 ±0.01 
19.28 ± 0.01 
20.62 ± 0.01 
18.86 ± 0.00 
20.56 ±0.01 
21.09 ±0.02 



18.90 ±0.01 



20.18 ±0.01 

21.01 ±0.01 
18.96 ±0.00 
20.98 ±0.01 
21.52 ±0.02 



2.1. HST versus ground-based photometry 

HST photometry for each cluster candidate is given in 
Table ED The HST magnitudes are given in the STMAG 
system, since transformation to the Johnson-Cousins sys- 
tem requires colour information which is not always avail- 
able. The photometry was obtained with the PHOT task 
in DAOPHOT, using an r = 0'.'5 aperture for the pho- 
tometry and an r = l'.'O — 1'.'5 annulus for the background 
measurements. The 0'/5 aperture contains about 90% of 
the light from a point source (e.g. Holtzman et al. [1995), 
but because the objects measured here are extended, an 
even larger fraction of the total light will fall outside the 
0'.'5 aperture. Thus, a total aperture correction of —0.2 
mag has been applied to the magnitudes in Table [SJ The 
exact correction will depend on the detailed spatial pro- 
hies of the clusters. In principle, this dependence might 
be estimated from the ISHAPE profile fits, but tying the 
photometry in Table |S] to the size measurements would 
make it difficult to backtrack these corrections. A con- 
stant —0.2 mag aperture correction is probably not far 
from the truth (Xarsen l2002[) . and more sophisticated cor- 
rections can easily be applied if desired. In cases where 



several observations were available for a cluster in a given 
band, Table|S]lists the measurement with the smallest for- 
mal error. 

In previous papers it has been documented that no 
systematic differences seem to be present between ground- 
based and HST-based colours for the cluster candidates 
(e.g. Larsen 120020 . although a random scatter of 0.1- 
0.2 mag exists. For integrated magnitudes, on the other 
hand, an offset of a few xO.l mag has been found be- 
tween ground-based and HST magnitudes, in the sense 
that ground-based data tends to give brighter magnitudes. 
With the larger sample of clusters available here, this com- 
parison can now be carried out in more detail. Of the clus- 
ters listed in Tabled 1245 have data in at least one of the 
filters F547M, F555W and F606W, all of which are rea- 
sonable approximations to the Johnson y-band. The mean 
difference (AVGround-WFPC2) between ground-based and 
WFPC2 photometry, including all objects with F547M, 
F555W or F606W HST data is -0.50 mag with a large 
scatter of er((AVGround-WFPC2}) = 0.74 mag. This scatter 
is partly due to the fact that the ground-based photom- 
etry of some of the fainter clusters has large errors, but 
decreases only slightly (to <r((AVGround-WFPC2)) = 0.65 
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mag) if clusters with formal errors larger than 0.2 mag on 
the ground-based V magnitudes are excluded. Thus, most 
of the errors are clearly of a systematic nature. 

If the samples observed with the Danish 1.54 m tele- 
scope and the NOT are compared with the HST photom- 
etry separately, interesting differences emerge. The im- 
age quality of the data taken with the two telescopes 
differ significantly, with typical FWHM seeing values of 
l'/5 and 0'.'8, respectively (Larsen 1999). For the galax- 
ies observed with the Danish 1.54 m, the mean dif- 
ference between ground-based and HST photometry is 
(AVcjround-WFPC2) = —0.76 ± 0.61 mag (where the 0.61 
mag refer to the standard deviation around the mean, 
not the error on the mean value). For the NOT sam- 
ple, the corresponding numbers are (AVGround-WFPC2) = 
—0.29 ± 0.50 mag. Thus, while the scatter remains large 
even in the NOT data, the systematic difference relative 
to the HST photometry is clearly smaller than for the 
Danish 1.54 m data. A few individual, relatively isolated 
clusters observed with both the NOT and Danish 1.54 m 
telescope have been analysed in detail (Larsen et al. 120011 
Larsen & Richtler 2004, in prep.) and for these clusters 
there is good agreement between ground-based and HST 
magnitudes (within ~ 0.1 mag). The differences between 
the mean magnitudes of the HST and ground-based sam- 
ples are therefore not due to trivial zero-point errors in 
the photometric calibrations. 

Clusters for which any of the comment flags in Tabled 
are set might be expected to show poorer agreement with 
the ground-based data. Indeed, if such clusters are rejected 
then (AVGround-WFPC2) = —0.40 ± 0.61 mag for all clus- 
ters, and (AVGround-WFPC2) = —0.59 ± 0.59 mag and 
(AVGround-WFPC2) = —0.10 ± 0.42 mag for the Danish 
1.54 m and NOT samples, respectively. Thus, the system- 
atic difference between HST and ground-based photome- 
try clearly decreases, albeit still with significant scatter. 
The remaining offsets can probably be attributed to con- 
tamination within the ground-based apertures which did 
not trigger any comment flags. In fact, if the HST pho- 
tometry is instead carried out using an r = l'.'O aperture 
radius, assuming that such an aperture encircles all flux 
from the objects (i.e. applying no aperture corrections) 
then the mean offset with respect to the NOT data is 
only (AVGround-WFPC2) = —0.03 ± 0.35 mag. Excluding 
clusters with comment codes, the difference decreases even 
further to (AVGround-WFPC2) = —0.003 ±0.33 mag. 

A smaller subset of the clusters have observations 
in HST bandpasses that allow a comparison with the 
ground-based colours. For example, Holtzman et al. 
(1995) give transformations to Johnson B — V colours 
for the F439W and F555W bandpasses, which are avail- 
able for 190 of the clusters in Table EJ The mean off- 
set between ground-based and HST B — V colours is 
(A(B — F)Ground-WFPC2) = —0.003 mag with a scatter 
of 0.17 mag. For the Danish 1.54 m and NOT samples, 
the differences are (A(B — F)Ground-WFPC2) = 0.05±0.15 
mag and (A(B-V) Ground- WFPC2) = -0.05 ±0.15 mag, 
respectively. This confirms that the ground-based colours 



are more accurate than magnitudes, presumably because 
the objects that contaminate the ground-based apertures 
tend to have similar ages and colours to the clusters them- 
selves. 

2.2. Tests of the profile-fitting algorithm 

Tests of the ISHAPE code have been carried out in several 
previous papers, in particular in Paper II. The reader is 
referred to that paper and to the documentation included 
with the code for further details. However, because the 
ability to fit the EFF a parameter is a more recent addi- 
tion, some tests of this particular feature are presented in 
the following. 

First, artificially generated sources with known pro- 
files were added to images of NGC 7793 and NGC 5194. 
The frames used for these tests were the WF4 chips from 
proposal 8591 (NGC 7793, F547M, 4x400 s) and proposal 
7375 (NGC 5194, F555W, 2x600 s). The artificial sources 
were generated by convolving EFF profiles with the HST 
PSF, generated by TinyTim, and then added to the sci- 
ence images with the MKSYNTH program (Paper II) . For 
each galaxy, artificial sources with FWHM of 1 and 2 pc, 
a = 1.0 and 1.5, and magnitudes of V — 20,21 and 22 
were added. For each combination of these input parame- 
ters, 25 objects were added at random positions and then 
fitted with ISHAPE using three different values (ao = 1.0, 
1.5 and 2.0) for the initial guess of a. 

The test results are summarised in TableEJ For brevity, 
results are shown for a — 1.0 only, but the conclusions 
reached for a = 1.5 are essentially identical. The first four 
columns in each row give the input V magnitude of the 
test objects, the FWHM in pixels and pc, and the a pa- 
rameter. The remaining columns list the mean S/N of the 
test objects (within the fitting radius), followed by the 
fitted FWHM and a values for the three initial guesses, 
ao, of the a parameter. The numbers in parantheses de- 
note the r.m.s. deviation around the mean, excluding the 
two most deviating data points at each extreme. The last 
two columns list the mean standard deviation ("03") of 
FWHM and a for the three fits to a single object. 

Neither the FWHM nor the a values returned by 
ISHAPE show any systematic dependence on ao, and the 
mean fitted FWHM and a are generally quite close to 
the input values. For S/N«50, the r.m.s. scatter of the 
FWHM and a measurements are typically 0.15-0.20 pix- 
els and 0.10-0.15 (dimensionless units), respectively. The 
mean standard deviation (03) of the three individual fits 
with different ao is much smaller than the object-to-object 
r.m.s. scatter, indicating that the uncertainty due to a par- 
ticular initial guess is smaller than the random measure- 
ment error. For typical cluster-like objects with S/N«50, 
these tests suggest that the FWHM and a parameter can 
be measured with an accuracy ~ 20% and 10-15%, re- 
spectively on a single image. 

Note, however, that the error on the effective radius 
R c ff derived from the FWHM and a can be very much 
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Table 6. Tests of the ISHAPE profile-fitting algorithm. For each combination of input parameters (V — 20/21/22, 
FWHM=lpc /2 pc) the output fitted FWHM and a parameters are shown for three initial guesses of a (ao = 
1.0/1.5/2.0). Numbers in parantheses denote the object-to-object r.m.s. deviation around the mean values, excluding 
the two most deviating points at each extreme. 





Input parameters 










Output parameters 
















a = 


1.0 




= 1.5 




= 2.0 






V 


FWHM, pix (pc) 


a 


(S/N) 


(FWHM) 


(a) 


(FWHM) 


(<*) 


(FWHM) 


(a) 


FWHM 


a 


NGC 77!):! (pS.-,<)l. F547M) 






















20.0 


0.62 (1.0) 


1.0 


95.4 


0.65 (0.10) 


1.01 (0.06) 


0.67 (0.10) 


1.01 (0.05) 


0.64 (0.10) 


1.00 (0.06) 


0.039 


0.023 


20.0 


1.25 (2.0) 


1.0 


96.7 


1.27 (0.14) 


0.99 (0.08) 


1.29 (0.15) 


1.00 (0.08) 


1.27 (0.13) 


0.99 (0.07) 


0.033 


0.019 


21.0 


0.62 (1.0) 


1.0 


40.0 


0.79 (0.15) 


1.07 (0.10) 


0.76 (0.16) 


1.06 (0.11) 


0.80 (0.21) 


1.08 (0.13) 


0.076 


0.044 


21.0 


1.25 (2.0) 


1.0 


41.1 


1.25 (0.25) 


1.01 (0.13) 


1.25 (0.27) 


1.00 (0.13) 


1.22 (0.31) 


0.99 (0.15) 


0.086 


0.044 
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larger, especially if a « 1. It is therefore important that 
any error selection be carried out on FWHM and a, and 
not on the derived quantity i? ff . 

As a final caution, the tests carried out here do not 
take into account any uncertainty on the PSF itself. In 
practice, this can be important even for HST images, since 
the so-called "diffusion" kernel makes an important con- 
tribution to the scattering of light between neighbouring 
pixels. While this effect is included in the modelling done 
by ISHAPE, the diffusion kernel has so far only been prop- 
erly characterised for the F555W filter. It is probably ap- 
propriate also for the F547M and F606W filters, but for 
bluer or redder passbands the true diffusion kernel could 
be significantly different. 

3. Contamination in the ground-based sample 

One of the motivations for this study was to test the relia- 
bility of the ground-based cluster identifications and quan- 
tify how much any mis-identifications might affect the spe- 
cific luminosities for the cluster systems derived in Paper 
III. To this aim, Table[7|lists photometry and object clas- 
sifications for all cluster candidates originally detected in 
the ground-based surveys, which are also included in the 
HST datasets. Of the 330 objects listed in Table □ 192 
are also included in Table OH The remaining 138 objects 
were excluded from the second round of fits (and Table|3J 
because of too low S/N for the variable-a fits. However, 
size information is listed for all objects in TableEJbased on 
the first round of profile fits (with a — 1.5). Because only 
one fit was done for each exposure, meaningful error esti- 
mates are only available for cluster candidates with 2 or 
more detections. The mean and median error on the clus- 
ter FWHM values for objects with multiple measurements 
are 0.5 pc and 0.26 pc, respectively, or 0.1-0.2 pixels (for 
typical galaxy distances of 5 Mpc) . This is consistent with 
the estimates of the accuracy of ISHAPE fits in Sec. 12.21 
(see also Paper II). 

The distribution of FWHM values for objects listed in 
Table is shown in Fig. 03 Only 6 objects were classified 
as 'likely/certain non-clusters' during the visual inspec- 



Table 8. Contamination of ground-based sample for indi- 
vidual galaxies. N to t is the total number of cluster candi- 
dates identified in the ground-based surveys, Nhst is the 
subset of those candidates covered by HST images and 
N<o.5 pc is the number of unresolved sources. The adopted 
distance moduli are also listed (see Paper I for references) . 
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0.86 


NGC 1313 


28.2 


45 


20 


2 


0.81 


0.68 


NGC 2835 


28.9 


9 


8 


6 


0.20 


0.17 


NGC 2997 


29.9 


34 


8 


1 


0.96 


0.98 


NGC 3184 


29.5 


13 


9 


3 


0.79 


0.83 


NGC 3621 


29.1 


45 


29 


6 


0.79 


0.75 


NGC 4258 


29.47 


44 


26 


2 


0.91 


0.91 


NGC 5055 


29.2 


24 


13 


3 


0.77 


0.76 


NGC 5194 
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1.00 
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30.5 
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1.00 


1.00 


NGC 7793 


27.6 


20 


11 
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0.46 
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tion, but 69 objects (or about 21% of the sample) have 
FWHM<0.5 pc and are thus essentially unresolved even 
on the HST images. Unresolved objects are not necessarily 
individual, isolated stars, but can also be loose groupings 
of stars (OB associations etc.) where ISHAPE simply picks 
one star and fits it. This is the case for about 25 out of 
the 69 unresolved objects, or about 1/3. A few examples 
are shown in Fig. [7| The size distribution in Fig. \E\ does 
not appear strongly peaked at a particular value (except 
for the unresolved sources near FWHM=0), but spans a 
range from the resolution limit up to 8-10 pc. If objects 
with FWHM<0.5 pixels are excluded, the formal estimate 
of the mean FWHM is 3.8 pc. For an EFF model with 
a = 1.5, this corresponds to a half-light radius of 4.3 pc (or 
4.0 pc if the profile is truncated at 50 pc), consistent with 
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Table 7. Clusters from ground-based survey reidentified in HST images. Photometry is from ground-based data. 
FWHM is the full-width-at-half-maximum in pc derived from EFF a — 1.5 fits. The first 10 rows of the table are 
reproduced here; the full table (330 rows) is only available in electronic form. 
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Fig. 5. Size distribution for cluster candidates identified 
in ground-based surveys. 



the typical sizes of young star clusters in the Milky Way 
and elsewhere. The mean FWHM changes only slightly 
(to 3.9 pc) if objects with 'Fit OK? = N' are excluded. 
The issue of cluster sizes will be discussed in more detail 
below for the full sample. 

Figure HO shows the distributions of absolute My 
magnitudes for unresolved and resolved objects. While 
the magnitude distribution for resolved objects (bottom 
panel) does extend to brighter magnitudes than for unre- 
solved ones, there are several unresolved objects brighter 
than My = —10. Such bright objects are unlikely to be 
individual stars, but most of them have the 'c' comment 
set (companions within l'/5), indicating that the ground- 
based magnitudes are likely contaminated by nearby ob- 
jects. Another possibility is that some of them are very 
compact star clusters. At magnitudes fainter than My = 
—9, Fig. shows a clear excess of unresolved objects (26 
out of 66, or 39%, compared to 21% for the whole sample), 
many of which may indeed be individual stars. 

Another way of estimating the contamination fraction 
is to use the object types from the visual inspection. Of the 
330 objects, 61 are labelled as type '2' (uncertain) or '3' 
(very likely non-cluster). Although classification as type 
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Fig. 6. Magnitude distributions for ground-selected clus- 
ter candidates with FWHM<0.5 pc (top) and FWHM>0.5 
pc (bottom). Sizes are from EFF fits with fixed a = 1.5. 
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Fig. 7. Examples of objects identified as clusters from 
the ground, but which are unresolved according to the 
ISHAPE fits 



'2' does not necessarily mean that the object in question 
is a non-cluster, these numbers again suggest a contami- 
nation rate of order 20% or less. 
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Fig. 8. Specific ?7-band luminosity (Tl(U)) vs. area- 
normalised star formation rate (Ssfr) for original ground- 
based data and updated values (filled diamonds) based on 
contamination rates estimated from HST data (Table |HJ- 

How does contamination of the cluster samples affect 
the relations between specific luminosity (Tl(U)) and host 
galaxy properties derived in Paper III? If the contamina- 
tion rate were the same in different galaxies, all specific 
luminosities would just be scaled down by a constant fac- 
tor, leaving any relations involving T^{U) intact. However, 
it is possible that the ratio of clusters to potential contam- 
inants varies from galaxy to galaxy. Table |H1 lists the total 
number of cluster candidates identified from the ground 
(Ntot) in each galaxy, followed by the number of ground- 
identified candidates contained within the WFPC2 images 
covering that galaxy (Nhst), and the number of objects 
among these that are unresolved (N<o.5 pc )- Based on the 
ground-based photometry, the V- and [/-band luminosi- 
ties in resolved objects, relative to all objects in the HST 
images, are listed in the last two columns. With a couple 
of exceptions, exclusion of unresolved objects affects the 
luminosities of the cluster systems only weakly. The two 
most notable exceptions are NGC 2835 and NGC 5204. In 
Paper I, we noted that all 7 clusters in NGC 5204 are quite 
red and resulted in an unusually high specific frequency 
(number of clusters relative to the host galaxy luminos- 
ity). Now, it appears that most of these objects may in 
fact be foreground stars. 

An updated version of the Tl(U) vs. Esfr plot from 
Paper III is shown in Figure |H1 The original data are 
shown with error bars, while new updated Tl(U) val- 
ues, corrected according to Table [SJ are shown with filled 
diamonds. The basic trend for Tl(U) to increase as a 
function of Ssfr is clearly preserved. The scatter in- 
creases somewhat after the exclusion of unresolved ob- 
jects. However, with only a couple of clusters left in galax- 
ies like NGC 2835 and NGC 5204 after exclusion of unre- 
solved objects, Tl(U) as defined in Paper III is probably 
no longer a good approximation to the true relative lu- 
minosity of the cluster system. To obtain a more useful 
number, one would likely have to probe to fainter mag- 



nitudes and thereby sample the cluster population more 
completely. 

4. The full sample 

4.1. Cluster ages 

The broad-band colours of simple stellar populations (such 
as star clusters) are functions of both age and metallicity, 
with additional complications arising from unknown red- 
denings and stochastic effects due to the finite number 
of stars in a cluster (e.g. Girardi et al. I19B5J1 . However, 
for clusters younger than ~ 10 9 years it is still possi- 
ble to obtain reasonably accurate photometric age esti- 
mates, especially if [/-band data are included. For such 
clusters, metallicity effects are weak, except for a brief pe- 
riod around 10 7 years when the integrated light is domi- 
nated by red supergiants. Here, cluster ages were obtained 
by fitting Bruzual & Chariot (2000; priv. comm.) model 
colours to the observed UBVI cluster colours. The SSP 
model fits were done by minimizing the r.m.s. deviation 
between model- and observed colours (weighted by their 
errors) as a function of age and reddening. In order to 
reduce the uncertainty on the age determinations, age es- 
timates were made only for clusters with <jb-v < 0.2 mag, 
(Tv-i < 0.2 mag and <tu-b < 0.3 mag. 

The ages of individual clusters derived from broad- 
band colours should only be regarded as approximate. The 
ground-based apertures may be contaminated by objects 
other than the cluster candidate itself, and model uncer- 
tainties also make the absolute ages uncertain. However, 
the relative age ranking of clusters should still be rea- 
sonably reliable. Line emission can also affect the broad- 
band colours of very young objects, and must be taken 
into account if accurate age estimates for objects younger 
than ~ 10 7 years are required (e.g. Anders & Fritze-v. 
Alvensleben l2UU3|) . 

4.2. Cluster sizes 

Table lists the mean FWHM and R c g for clusters in 
each galaxy. Only clusters classified as 'Type 1' and with 
'Fit OK = YES' were included. The number of clusters 
in each galaxy satisfying these criteria is given in the sec- 
ond column. Because of the poorly defined R c h for clus- 
ters with a < 1, such clusters were excluded before com- 
puting the average i? c ff values in Table Objects with 
very steep envelope slopes (a > 5) represent clear out- 
liers (see Section fOl below) and were also excluded from 
further analysis. The mean values given in the table are 
unweighted. This is because there is a correlation between 
the cluster sizes and their associated errors (the relative 
errors remain roughly constant), which would lead to a 
strong bias in the mean if the size measurements were 
error- weighted . 

The mean effective radius is (R e s) — 3.94 ± 0.12 pc, 
perhaps slightly larger compared to those for the ground- 
selected sample (SecOJ and other young and old star clus- 
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Fig. 9. Mean FWHM (top) and effective radii (bottom) 
for stellar clusters versus parent galaxy distance modulus. 
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Fig. 10. Mean FWHM (top) and effective radii (bottom) 
for stellar clusters versus parent galaxy Ssfr (derived 
from IRAS far-infrared luminosities, as described in Paper 
III). 

Table 9. Mean FWHM and effective radii for clusters in 
galaxies. 

Galaxy N (FWHM) (R eB ) (a > 1) (R cff ) (all) 

pc pc pc 



ters. This may be due to the fact that the clusters are as- 
sumed to follow a single power-law out to a total radius of 
50 pc, while in reality the behaviour at large radii is poorly 
constrained. If the cluster profiles decline more rapidly at 
large radii the effective radii would decrease, especially 
for objects with a- values close to 1. Also, the size cut im- 
posed in order to exclude point sources may introduce a 
bias against the most compact clusters. In the last column 
of Table E| the a > 1 requirement is abandoned. Clearly, 
this leads to an increase in the mean effective radii, but it 
is stressed that the R c s values for these objects are very 
uncertain and depend strongly on the adopted outer radii. 

The mean FWHM and R cS from Table H are plot- 
ted versus parent galaxy distance modulus in Fig. 03 The 
sizes do show some correlation with parent galaxy dis- 
tance, possibly due to a less than perfect correction for 
the PSF. It is also possible that a larger fraction of the 
objects detected in more nearby galaxies are individual 
stars, rather than clusters, which made it into the list of 
cluster candidates despite the size cut. Contributing to 
this effect, the number of individual stars bright enough 
to be detected in the ground-based photometry would in- 
crease at smaller distances. Furthermore, the fixed size cut 
at FWHM=0.2 pixels corresponds to a different physical 
cluster size in different galaxies, ranging from 0.3 pc (core 
radius ~ 0.15 pc) in NGC 7793 to 1.2 pc (core radius ~ 0.6 
pc) in NGC 7424. 

Fig. EH shows the mean cluster sizes as a function of 
the area-normalised host galaxy star formation rate, Esfr 
(see Paper III). The scatter in the FWHM plot is some- 
what larger than for the i? e ff values, but neither FWHM 
nor i? c ff shows any obvious correlation with Xsfr- Below, 
the data for all galaxies is combined and analysed collec- 
tively in order to improve statistics, but in order to re- 
duce possible systematic effects due to different distances, 
clusters in the closest (NGC 7793) and two most distant 
galaxies (NGC 2997 and NGC 7424) are excluded for the 
analysis of structural parameters. Thus, in summary the 
selection parameters applied to the cluster candidates in 
Table for the following analysis are: 

— Object type = 1 (certain / very likely cluster) 

— Fit = OK 

— When i? e ff is involved: require 1 < a < 5 

— Photometric errors: (Jb—v < 0.2 mag, crv-i < 0.2 mag 
and ajj-B < 0.3 mag 

— Exclude NGC 2997, NGC 7424 and NGC 7793 

Figure ^2 shows the distributions of FWHM and i? c ff 
values for the combined sample. The paucity of objects 
with FWHM< 1 pc is partly a selection effect, due to 
the size cut at FWHM=0.2 pixels, so it is not clear that 
the distribution of FWHM even has a well-defined peak. 
However, most clusters have FWHM<6-8 pc, similar to 
the range spanned by the ground-selected clusters (Fig.|SJ. 
The i? c ff distribution is more complicated to interpret be- 
cause the cut in FWHM does not translate to a unique 
i? e ff but depends on a. For a typical a ~ 1.3 (see below), 
a FWHM of 1 pc corresponds to R^g = 1.6 pc, but this 
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Fig. 11. Distributions of FWHM and effective radii for the 
combined cluster samples. In the lower panel, the hatched 
and outlined histograms are for clusters with 1 < a < 5 
and a < 5, respectively. 



increases to i? c ff = 3 pc for a = 1.1 and i? e ff = 8 pc 
for a = 0.9. However, there is a fairly rapid drop-off in 
the distribution for i? c g > 5 pc, and the vast majority 
of the clusters have i? e ff < 10 pc. The mean FWHM and 
i? e ff values given in Table El and shown in Figs. l9l and ITUI 
should clearly be interpreted with caution. 

4.3. Envelope slope versus age 

Figure IT21 shows the distribution of a values for clusters 
in four different age bins. For reference, typical errors of 
0.15 on a are indicated (cf. Sec 12. 2(1 . The distribution of a 
values is peaked around 1.2—1.3 in all four bins, but there 
may be a trend of mean a increasing with age, at least 
in the sense that the youngest age bin appears to contain 
a higher fraction of clusters with a < 1. One important 
caveat in this comparison is that the mean mass is likely 
to increase from the youngest to the older age bins, due 
to the increase in M/L ratio with age. Thus, in principle 
the difference between the a-distributions in Fig. ll2l might 
be due to the different mass ranges sampled in each bin, 
rather than being an evolutionary effect. To test whether 
this might be the case, Fig.^Jdisplays the a-distributions 
in two age intervals, but now also divided into different 
mass intervals. The number of clusters in each panel is 
quite small, but the excess of clusters with a < 1 seems 
to be present for log(age)<7.0 in all three mass bins. This 
suggests that the difference between the a-distributions in 
Fig- El is n °t j us t an effect of the different mass intervals 
covered at different ages. 
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Fig. 12. Distribution of envelope slopes in four age bins 

Other authors have previously noted that very young 
clusters tend to be surrounded by relatively extended en- 
velopes with more light at large radii than a King profile. 
Elson et al. (1987) estimated that as much as 50% of the 
mass in young LMC clusters may reside in unbound halos. 
In the Milky Way, a large fraction of the youngest open 
clusters also have very large radii and may be unbound 
and in the process of dispersing away (Janes et al. [1988). 
Whitmorc et al. (1999| showed radial profiles for three 
clusters in the Antennae, illustrating a gradual transition 
from extended envelopes with no well-defined outer radius 
(for the highly luminous 'Knot S', only a few Myr old) to 
older clusters where a tidal cut-off becomes apparent. A 
similar extended halo was observed for a very luminous, 
15 Myr old cluster in NGC 6946 (Larsen et al. lSUUT)) . The 
structure of very young clusters may hold important clues 
to the structure of the progenitor clouds out of which the 
clusters formed, although it may prove challenging to dis- 
entangle this from the effects of early dynamical evolution. 

4.4. Cluster sizes versus age and mass 

Figures ITU and IT51 show FWHM and i? c ff versus age and 
mass for the combined cluster samples. Clusters older than 
1 Gyr are excluded from the plots because of large un- 
certainties on the ages (and therefore also on the masses 
derived from integrated photometry). The masses were es- 
timated using M/L ratios from the Bruzual & Chariot 
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Fig. 13. Envelope slopes in different mass bins for two 
different ages 
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Fig. 14. Cluster FWHM versus age (top) and mass (bot- 
tom). The dashed curve in the lower panel represents a 
relation corresponding FWHMocM 1 / 3 (constant density). 
The solid line is a least-squares power-law fit to FWHM 
vs. mass. 
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Fig. 15. Effective radius (R c g) versus age (top) and mass 
(bottom panel) . Clusters with a < 1 and a > 1 are shown 
with diamonds and plus markers, respectively. The dashed 
curve in the lower panel represents a relation correspond- 
ing i? c ff ocM 1 / 3 (constant density). The solid line is a least- 
squares power-law fit to R c s vs. mass for clusters with 
a > 1. 



models, assuming a Salpeter IMF from 0.1 Mq— 100 Mq. 
Absolute masses are sensitive to the shape and lower-mass 
cut-off of the IMF, and would be ~ 30% lower for a Kroupa 
(2002) IMF, but this is not important for the relative com- 
parison attempted here. In Fig. ^1 clusters with a < 1 
and a > 1 are shown with different symbols (diamonds 
and plus markers). The dashed lines superimposed on the 
lower panels (size vs. mass) of each figure represent the 
size oc M 1 / 3 relation corresponding to constant cluster 
density. Typical error bars are also shown in the top right 
corner in each of the lower plots. 

Neither Fig. 1141 nor Fig. 1151 shows any strong evidence 
for a general correlation between cluster size and age. 
Since FWHM is roughly equivalent to 2 x r c , it is interest- 
ing to note that Mackey & Gilmore ( 2003 ) found a strong 
correlation between core radius and age for clusters in the 
LMC (confirming earlier results by Elson et al. 1989). For 
ages < 10 8 years, there are essentially no clusters in the 
LMC with r c > 2.5 pc, while clusters with ages ~ 10 9 
years show a full range of core radii from less than 1 to 8 
pc. While it cannot be ruled out that the most extended 
clusters are missing from the sample, Fig. 1141 shows no 
strong r c -age relation similar to that in the LMC. 

Likewise, the size-mass plots show no strong correla- 
tions. The solid line in the bottom panel of Fig. [21 is a 
least-squares power-law fit of the form 



FWHM = A(M/MqY 



(2) 
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Fig. 16. Histograms of the R c s distributions in four mass 
ranges. Filled histograms are for clusters with a > 1, while 
the outlined histograms include clusters with a < 1. The 
shaded areas indicate the range of sizes in each mass in- 
terval corresponding to an i? e ff oc M 1 / 3 relation. 

where M is the mass and FWHM is measured in pc. The 
best-fitting parameter values are A = 0.94 ± 0.25 pc and 
B = 0.08 ±0.03. The fit formally indicates a ~ 3er correla- 
tion between FWHM and mass, but with a large scatter. 
The errors on the coefficients A and B are the formal er- 
rors from the least-squares fit, but a jack-knife test yields 
nearly identical error estimates. The size cut may intro- 
duce a bias against more compact, low-mass objects, but 
the fit does not change substantially even if restricted to 
objects more massive than 10 4 M Q . While no strong con- 
straints can be put on a FWHM vs. mass relation based on 
Fig. 1141 the formal relation derived from the least-squares 
fit is shallower than the M 1 / 3 relation corresponding to 
constant core density (indicated by the dashed line). A 
similar fit is carried out for the R c s vs. M data in Fig. 1151 

i? off = A(M/M G ) B , (3) 

but only for clusters with 1 < a < 5. Here, the best- 
fitting parameters are A — 1.12 ± 0.35 pc and B = 0.10 ± 
0.03, quite similar to the slope of the FWHM vs. mass 
relation in Fig. 1141 and again only significant at about the 
3cr level. Formally, the slope is slightly steeper than (but 
compatible with) the r oc L 07 relation found by Zepf 
et al. fiMty for young clusters in NGC 3256. As for the 



mean R e g values in Table an important caveat is that 
the derived i? e ff are sensitive to the somewhat arbitrary 
truncation of the luminosity profiles at 50 pc. The slope 
of the i? G ff vs. M relation remains stable at least for cut- 
offs between 30 and 100 pc, even if the clusters with a < 
1 are included in the fit (but the normalisation changes 
somewhat). As in the fit to FWHM vs. mass (and Table|5J), 
the individual datapoints are not weighted by their errors. 

Using the HST photometry in F547M, F555W or 
F606W in Table |5] instead of the ground-based data to 
estimate the U-band luminosities of the clusters has no 
significant effect on the FWHM vs. M and R e g vs. M 
relations. For i? e ff vs. M, the slope B decreases by only 
0.003, while for FWHM vs. M it decreases by 0.014, in 
both cases well within the formal uncertainties on the fits. 
Note, however, that the ages still have to be estimated 
from the ground-based photometry. 

The high degree of crowding in Fig. ^] makes it diffi- 
cult to visually assess to what extent the two curves agree 
with the data. A somewhat different representation of the 
same data is shown in Fig. 1161 where the R c g distributions 
are shown for five separate mass bins. The sample with 
a > 1 is shown with shaded histograms, while the sam- 
ple including clusters with a < 1 is shown with outlined 
histograms. In each panel, the range of sizes correspond- 
ing to an i? e ff oc M 1//3 scaling is also indicated (using the 
same arbitrary normalisation as in Fig. 1151 not taking any 
scatter into account). The paucity of objects in the 0-1 pc 
bin is due to the size cut, but even if one accounts for the 
fact that the most clusters may have been systematically 
excluded, it seems difficult to reconcile the overall i? c ff 
distributions in Fig. 1161 with a constant-density relation. 

Although completeness effects are difficult to quantify, 
they would most likely tend to work against the detection 
of faint, extended clusters, thereby strengthening rather 
than weakening any existing size-mass trend. Thus, while 
the size-mass trends suggested by Figs. 1141 and 1151 remain 
quantitatively uncertain, the current data seem to imply 
that any size-mass trend, if it exists, is substantially shal- 
lower than for a constant-density relation. 

4.5. Environment 

While more subjective than the measurements of FWHM, 
a and i? c ff , the comment flags in Table hold useful infor- 
mation about the surroundings of each cluster. It is reiter- 
ated that the comment flags are based on visual inspection 
of a fairly heterogeneous dataset, and that the morphology 
of cluster candidates can be quite wavelength-dependent. 
Note, however, that the comment codes are based on vi- 
sual inspection of all available images of each cluster, and 
the majority of the clusters have imaging in a B-, V- or 
J-band equivalent filter. Table ITU1 lists the fraction of clus- 
ter candidates with comment flags a/c/d (all of which are 
likely indicators of multiplicity) and e/f/g/h (more gen- 
eral crowding indicators). These statistics are given both 
for the best cluster candidates with structural parameters 
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Table 10. Statistics on comment flags in different age intervals 



Age range 


Fit-OK 


= YES 




All 






Type 


= 1 










a/c/d 


e/f/g/h 


Type=2 


a/c/d 


e/f/g/h 


6.0 < log(age) < 7.0 


36% 


2% 


42% 


54% 


19% 


7.0 < log(age) < 7.5 


22% 


0% 


22% 


32% 


4% 


7.5 < log(age) < 8.0 


16% 


0% 


25% 


29% 


7% 


8.0 < log(age) < 9.0 


10% 


0% 


7% 


14% 


0% 



determinations (i.e. Fit-OK = YES and Type 1), as well 
as for all potential cluster candidates including those with- 
out reliable fits (Fit-OK = YES/NO and Type 1/2). Data 
for all galaxies are included in Table but the num- 
bers remain unchanged within 1-2% even if NGC 2997, 
NGC 7793 and NGC 7424 are omitted as in the previous 
sections. 

Very few clusters with Fit-OK=YES have any of the 
/e/f/g/h flags set. This is no coincidence, because these 
flags indicate exactly those conditions which would make 
profile fits uncertain. The table indicates a strong evo- 
lution in the environment as a function of age. Most 
of the cluster candidates in crowded environments (flags 
e/f/g/h) are younger than 10 7 years. The tendency for the 
crowding to decrease with age is probably a consequence 
of fading and dispersion of the surrounding stellar popula- 
tion. Assuming typical velocity dispersions of a few km/s 
within (unbound) star forming regions, the expansion will 
amount to a few tens of pc in 10 7 years. 

Of the objects in the youngest age bin, 42% were clas- 
sified as "uncertain" (type 2) . This underscores the funda- 
mental difficulty of identifying the youngest clusters. The 
problem of defining an appropriate selection criterion for 
bona-fide clusters is far from trivial. Sometimes the main 
problem is simply that an object is only barely resolved. 
In such cases, better angular resolution would help con- 
firm or rule out the cluster nature. For objects with com- 
plex morphology it can be difficult to determine whether 
or not an object is a true star cluster, even if resolution 
would otherwise not be a problem. Examples can be seen 
in Fig. 01 and in panels (e), (f) and (g) of Fig. [21 In these 
cases, it is difficult to determine whether a well-defined 
stellar cluster is present. For low-mass clusters of low age, 
an additional problem is that the integrated light can be 
dominated by a few luminous stars, making it difficult to 
distinguish such objects from random superpositions of 
individual stars along the line-of-sight. 

As pointed out in the introduction, most stars prob- 
ably form in clusters, but only a small fraction of young 
embedded clusters survive as bound objects. It may also 
happen that only a fraction of the stars in a cluster remain 
bound, while the rest disperse away (Kroupa 2001} . Thus, 
a few Myr-old concentration of stars may be a bound star 
cluster, a bound core surrounded by an expanding enve- 
lope, or an entirely unbound association which will soon 
disperse away completely. Other star formation may also 
be taking place nearby, perhaps triggered by the young 



cluster. So it is not surprising that a large fraction of the 
youngest objects have a messy morphology. 

The age distribution of double or multiple objects does 
not appear to be as strongly peaked at young ages, with 
some objects flagged 'a/c/d' even in the oldest age bin. It 
is possible that at least a fraction of these objects are gen- 
uine double clusters, similar to those found in the Large 
Magellanic Cloud. The LMC binary clusters tend to be 
predominantly young, though a few pairs as old as several 
times 10 8 years exist (Dieball et al. 2002). This seems to be 
consistent with the decreasing fraction of multiple objects 
at high ages in Table^U However, because the main source 
of photometry in this paper is ground-based imaging, no 
information is available about possible colour/age differ- 
ences between the components in such pairs. Multi-colour 
HST imaging, especially including [/-band data, would al- 
low a more thorough investigation of double clusters and 
make a comparison with the LMC sample possible. 

5. Summary and conclusions 

Using a combination of HST/WFPC2 imaging and 
ground-based UBVI photometry, a catalogue of cluster 
candidates in 18 nearby spiral galaxies has been compiled. 
Only objects with a S/N>50 on the HST images (within 
an r = 5 pixels aperture) were included, allowing for a 
detailed analysis of the structure of individual clusters. 
Analytic profile fits of the form P{r) cx (1 + (r/r c ) 2 )~ a 
were carried out, including a proper modelling of the 
HST/WFPC2 PSF, and allowing both the core radius 
r c and envelope slope parameter a to vary Each cluster 
candidate has also been visually inspected and comment 
flags relating to crowding and multiplicity are given. These 
comment flags, combined with the photometric data and 
structural parameters, may be helpful when using the list 
of cluster candidates to select targets e.g. for spectroscopic 
studies. 

The HST imaging indicates a mean contamination rate 
of 20% or less for the ground-based cluster surveys in 
Paper I— III, although the contamination rates in some in- 
dividual galaxies (most notably NGC 5204) are higher. 
However, the relation between specific [/-band luminos- 
ity of the cluster systems Tl(U) and Ssfr (Paper III) 
remains valid after correction for contamination. Because 
very few clusters remain after the contamination correc- 
tion in some galaxies, sampling the cluster populations to 
fainter magnitudes than the limits defined in Papers I— III 
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would probably reduce the scatter in the Tl{U) vs. Esfr 
relation. 

The cluster catalogue has been used to investigate 
trends and relations between various cluster properties, 
although the analysis is complicated by the fact that se- 
lection effects are difficult or impossible to control with 
a dataset as heterogeneous as the one used here. Most 
clusters have FWHM less than about 8 pc with a formal 
mean of about 2.7 pc, corresponding to a mean core radius 
(r c ) « 1.3 pc, but very compact clusters may be missing 
because of the size cut imposed in order to exclude point 
sources (individual stars). For the subset of the clusters 
which have 1 < a < 5 and therefore reasonably reliable 
measurements of the effective radius i? e ffi the mean value 
is (R e g) = 3.94 ± 0.12 pc, but this mean value may again 
be affected by a selection bias against the most compact 
clusters. The effective radii are also sensitive to the poorly 
constrained behaviour of the luminosity profiles at large 
radii. Here, the profiles are truncated at 50 pc. The struc- 
tural parameters show little or no variation from galaxy 
to galaxy, especially when considering that the distances 
are not always known very accurately. In particular, the 
effective radii are uncorrelated with the host galaxy area- 
normalised star formation rate, and are also very similar 
to those of open and globular clusters in the Milky Way, 
globular clusters in early- type galaxies, and young clusters 
in merger galaxies and starbursts. It is quite remarkable 
that the sizes of stellar clusters are largely invariant with 
respect to the properties of the parent galaxy. Physical 
parameters such as gas density and -pressure probably 
play a major role in regulating the star formation rate 
(Kennicutt 11998(1 . but while these factors may affect the 
formation efficiency of bound clusters (Paper III) they do 
not appear to have a strong impact on the structure of the 
clusters themselves, once formed. Exceptions to this rule 
do exist, including the "faint fuzzy" star clusters observed 
in some lenticular galaxies (Brodie & Larsen 2002), and 
there is a general trend for the sizes of globular clusters to 
increase as a function of galactocentric distance (van den 
Bergh et al. HMTjl 

While both the FWHM and i? e ff are found to corre- 
late with cluster mass, least-squares power-law fits yield 
slopes that are substantially shallower than for a constant- 
density relation, implying an increase in cluster density as 
a function of mass. Qualitatively, these results are in agree- 
ment with previous data for young star clusters as well as 
old globular clusters. Quantitatively, the relations show a 
large scatter and remain uncertain. Ashman & Zepf |2001 ) 
have argued that an increasing star formation efficiency as 
a function of cluster mass may explain (at least partially) 
the lack of a strong size-mass relation. 

Many of the youngest clusters have extended, shallow 
outer envelopes. This tendency seems to be a general one, 
noted previously for a few isolated cases in the Antennae, 
NGC 6946, and for LMC clusters. The structure of these 
young objects may hold important clues to the early dy- 
namical evolution of clusters and the density distribution 
of the parent proto-cluster clouds. Older clusters gradu- 



ally evolve towards King-type profiles with a finite tidal 
radius. 

Finally, a strong correlation between cluster age and 
crowding is found, with most of the strongly crowded 
clusters having ages < 10 7 years. About 1/3-1/2 of 
these young objects are double or multiple sources, but 
the identification as bona-fide clusters is often uncertain 
even on WFPC2 images in these fairly nearby galaxies. 
Future multi-colour imaging with the Advanced Camera 
for Surveys on HST should help resolve many of the dif- 
ficulties encountered in this study. 
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